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A Computational Fluid Dynamics 
Study of Transitional Flows in 
Low-Pressure Turbines Under a 
Wide Range of Operating 
Conditions 

A transport equation for the intermittency factor is employed to predict the transitional 
flows in low-pressure turbines. The intermittent behavior of the transitional flows is taken 
into account and incorporated into computations by modifying the eddy viscosity, with 

the intermittency factor, y. Turbulent quantities are predicted by using M enter’s two- 
equation turbulence model (SST). The intermittency factor is obtained from a transport 
equation model which can produce both the experimentally observed streamwise varia- 
tion of intermittency and a realistic profile in the cross stream direction. The model had 
been previously validated against low-pressure turbine experiments with success. In this 
paper, the model is applied to predictions of three sets of recent low-pressure turbine 
experiments on the Pack B blade to further validate its predicting capabilities under 
various flow conditions. Comparisons of computational results with experimental data 
are provided. Overall, good agreement between the experimental data and computational 
results is obtained. The new model has been shown to have the capability of accurately 
predicting transitional flows under a wide range of low-pressure turbine conditions. 
[DOI: 10.1115/1.2218888] 


1 Introduction 

The process of transition from laminar to turbulent flow is a 
major unsolved problem in fluid dynamics and aerodynamics. One 
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area where the transition process plays an important role and is 
even more complicated due to the diverse flow conditions encoun- 
tered is the low-pressure turbine applications. Transitional flows 
in these applications are affected by several factors such as vary- 
ing pressure gradients, wide range of Reynolds number and 
freestream turbulence variations, flow separation, and unsteady 
wake-boundary layer interactions. Accurate simulation and pre- 
diction of transitional flows under these diverse conditions is key 
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to design of more efficient jet engines. 

In low-pressure turbine applications, flow over the blades is 
mostly turbulent at the high Reynolds number conditions encoun- 
tered at takeoff and the efficiency is at its design maximum. How- 
ever, at lower Reynolds number conditions which correspond to 
high altitudes and cruise speeds the boundary layers on the airfoil 
surface have a tendency to remain laminar; hence, the flow may 
separate on the suction surface of the turbine blades before it 
becomes turbulent. This laminar separation causes unpredicted 
losses, substantial drops in efficiency, and increase in fuel con- 
sumption [1-3]. 

In order to calculate the losses and heat transfer on various 
components of gas turbine engines, and to be able to improve 
component efficiencies and reduce losses through better designs, 
accurate prediction of development of transitional boundary layers 
is essential [1], 

One approach proven to be successful for modeling transitional 
flows is to incorporate the concept of intermittency into computa- 
tions. This can be done by multiplying the eddy viscosity obtained 
from a turbulence model, fi„ used in the diffusive parts of the 
mean flow equations, by the intermittency factor, y (Simon and 
Stephens [4]). This method can be easily incorporated into any 
Reynolds averaged Navier-Stokes solver. In this approach, the in- 
termittency factor, y, can be obtained from an empirical relation 
such as the correlation of Dhawan and Narasimha [5], or it can be 
obtained from a transport model. 

Dhawan and Narasimha [5] correlated the experimental data 
and proposed a generalized intermittency distribution function 
across flow transition. Gostelow et al. [6] extended this correlation 
to flows with pressure gradients under the effects of a range of 
ffeestream turbulence intensities. Solomon et al. [7], following the 
work of Chen and Thyson [8], developed an improved method to 
predict transitional flows involving changes in pressure gradients. 
These empirical methods led to development of transport equa- 
tions for intermittency. 

Steelant and Dick [9] proposed a transport equation for inter- 
mittency, in which the source term of the equation is developed 
such that the y distribution of Dhawan and Narasimha [5] across 
the transition region can be reproduced. Steelant and Dick used 
their model, coupled with two sets of conditioned Navier-Stokes 
equations, to predict transitional flows with zero, favorable, and 
adverse pressure gradients. However, since their technique in- 
volved the solution of two sets of strongly coupled equations, the 
method is not compatible with existing computational fluid dy- 
namics (CFD) codes, in which only one set of Navier-Stokes 
equations is involved. Moreover, the model was designed to pro- 
vide a realistic streamwise y behavior but with no consideration 
of the variation of y in the cross-stream direction. 

Cho and Chung [10] developed a k-e-y turbulence model for 
free shear flows. Their turbulence model explicitly incorporates 
the intermittency effect into the conventional k-e model equations 
by introducing an additional transport equation for y. They ap- 
plied this model to compute a plane jet, a round jet, a plane far 
wake, and a plane mixing layer with good agreement. Although 
this method was not designed to reproduce flow transition, it pro- 
vided a realistic profile of y in the cross-stream direction. 

Suzen and Huang [11] developed an intermittency transport 
equation combining the best properties of Steelant and Dick’s 
model and Cho and Chung’s model. The model reproduces the 
streamwise intermittency distribution of Dhawan and Narasimha 
[5] and also produces a realistic variation of intermittency in the 
cross-stream direction. This model has been validated against Eu- 
ropean Research Community On Flow Turbulence And Combus- 
tion (ERCOFTAC) benchmark T3-series experiments reported by 
Savill [12,13], low-pressure turbine experiments of Simon et al. 
[14], and separated and transitional boundary layer experiments of 
Hultgren and Volino [15] with success [11,16-21], 

In this paper we concentrate on prediction of three recent low- 
pressure turbine experiments on the Pratt and Whitney’s Pack B 
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Chord length, L 
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Axial chord to chord ratio, L/L=0.906 

Pitch to chord ratio, P/L = 0.8 

Blade inlet angle, /?, = 35° 

Blade outlet angle, = -60° 


Fig. 1 P&W Pack B blade cascade details 


blade under low Reynolds number conditions using the transport 
model for intermittency. Due to the fact that the Pack B blade is 
vety sensitive to changes of flow conditions, it is an ideal test 
blade for validating the transition/turbulence models. The three 
sets of experiments considered are conducted by Lake et al. 

[3.22] , Huang et al. [23], and Volino [24] at three independent 
facilities. These experiments provide an extensive database for 
investigating transitional flows under low-pressure turbine condi- 
tions and are employed as benchmark cases for further testing of 
the predicting capabilities of the current intermittency model. A 
summary of the experiments are given in the next section. In Sec. 
3, the intermittency transport model is presented and implementa- 
tion of the model and the empirical correlations employed for the 
onset of transition are described. In Sec. 4, the predictions of the 
new intermittency model are compared against the experimental 
data. Conclusions are provided in Sec. 5. 

2 Low-Pressure Turbine Experiments 

In this paper, we concentrate on computation of three sets of 
low-pressure turbine experiments using the intermittency transport 
model. These experiments are conducted by Lake et al. [3,22], 
Huang et al. [23], and Volino [24], In these experiments Pratt and 
Whitney’s Pack B blade is used; the details of the blade are shown 
in Fig. 1. Overall, these experiments cover a Reynolds number 
range from 10,000 to 172,000 and the ffeestream turbulence in- 
tensity ranges from 0.08% to 4%. The cases and data used for 
comparison in this paper are summarized in Table 1. In the fol- 
lowing sections details of these experimental efforts are given. 

2.1 Pack B Blade Cascade Experiments of Lake et al. 

[3.22] . Lake et al. [3,22] conducted experiments on the Pack B 
blade in order to identify methods for reducing separation losses 
on low-pressure turbine blades under low Reynolds number con- 
ditions. In the experiments, they investigated flows at low Rey- 
nolds numbers of 43,000, 86,000, and 172,000 based on inlet 
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Table 1 Details of the experiments used for comparison with computations 




C x 

Re 


Data used for 

Source 

Test Section 

(m) 

(UJCJv) 

FSTI (%) 

Comparison 

Lake et al. [3,22] 

P&W Pack B cascade 

0.1778 

86,000 

1 & 4 

C p distribution 




172,000 

1 & 4 

C p distribution 

Huang et al. T231 

P&W Pack B cascade 

0.1595 

10,000 

0.08 

C p distribution 




25,000 

0.08 

C p distribution 




50,000 

0.08, 1.6, 2.85 

C p distribution. 






velocity profiles 2 




75,000 

0.08, 1.6, 2.85 

C p distribution. 






velocity profiles 2 




100,000 

0.08, 1.6, 2.85 

C. distribution. 






velocity profiles 2 

Volino [24] 

P&W Pack B single 

0.1537 

10,291 

0.5 

C p distribution. 


passage 




velocity profiles 




20,581 

0.5 

C p distribution. 






velocity profiles 




41,162 

0.5 

C p distribution. 






velocity profiles 




82,324 

0.5 

C p distribution. 






velocity profiles 


“Velocity profiles are available for FSTI=0.08% and 2.85% from experiments. 


velocity and axial chord and freestream turbulence intensities 
(FSTI) of 1% and 4%. These conditions are similar to those en- 
countered at high-altitude, low-speed flight of reconnaissance un- 
manned aerial vehicles used by USAF. 

In Lake’s experiments, surface pressure coefficients, boundary 
layer velocity, and turbulence profiles, total pressure loss data 
were obtained at FSTI=1% and FSTI=4%. The test setup shown 
in Fig. 2 included eight blades with axial chord of 0.1778 m 


(7 in.), and blade spacing of 0.1575 m (6.2 in.). The blades were 
numbered 1 through 8 starting from the inside bend. Boundary 
layer measurements were taken on blade 5 and surface pressures 
were measured around blades 4 and 6. In this paper, the Pack B 
blade experiments with Reynolds numbers of 86,000 and 172,000 
and freestream turbulence intensities of 1 % and 4% are computed 
and comparison of pressure distributions between experiments 
and computations are performed. 



Fig. 2 Experimental setup used by Lake et al. [3,22] 
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Distance (m) 


Fig. 3 Comparison of computed and experimental decay of 
turbulence for experiments of Huang et al. [23], with grid 0 


2.2 Pack B Blade Cascade Experiments of Huang et al. 
[23]. Huang et al. [23] conducted experiments on Pack B blade 
cascade for a range of Reynolds numbers and turbulence intensi- 
ties. The Reynolds numbers range from 10,000 to 100,000 based 
on inlet velocity and axial chord as listed in Table 1. In their 
experiments the blades had an axial chord length of 0.1595 m 
(6.28 in.). The freestream turbulence intensity in the tunnel was 
measured as 0.08%. In order to increase the turbulence intensity, 
two grids with different mesh sizes were used. One of the grids 
had the mesh size of 0.0254 m (denoted as grid 0) and the other 
had 0.008 m (denoted as grid 3). The decay of turbulence after the 
grids was measured using crosswire and they are shown in Figs. 3 
and 4 along with the computed results for grid 0 and grid 3, 
respectively. The grids were movable in the tunnel so that the 
turbulence level of the flow that reaches the blades could be con- 
trolled by moving the grid that is, by increasing or decreasing the 
distance between the grid and the blade. Experiments were per- 
formed for Reynolds numbers 50,000, 75,000, and 100,000, with 
grids placed 0.762 m (30 in.) away from the blade leading edge, 
corresponding to turbulence intensities of 2.85% and 1.6% at the 
leading edge for grid 0 and grid 3, respectively. For Re 
= 100,000, grid 0 is placed at 0.5588 m (22 in.) and 0.3556 m 
(14 in.), corresponding to turbulence intensities of 3.62% and 
5.2%, respectively. Pressure coefficient data are available for all 
cases and detailed boundary layer measurements are available for 
Re =50,000, 75,000, and 100,000 with FSTI=0.08% and 2.85% 
cases. The cases and data used for comparisons in this paper are 
listed in Table 1. 



Distance (m) 


Fig. 4 Comparison of computed and experimental decay of 
turbulence for experiments of Huang et al. [23], with grid 3 



Fig. 5 Schematic of the test section for experiments of Volino 
[24] 


2.3 Pack B Experiments of Volino [24], Volino [24] investi- 
gated the boundary layer separation, transition, and reattachment 
under low-pressure turbine airfoil conditions. The experiments in- 
cluded five different Reynolds numbers ranging between 10,291 
and 123,492 and freestream turbulence intensities of 0.5% and 
9%. The test section consisted of a single passage between two 
Pack B blades as shown in Fig. 5. The axial chord length of the 
blades was 0.1537 m (6.05 in.). There are flaps located upstream 
of each blade to control the amount of bleed air allowed to escape 
from the passage. These flaps were adjusted by matching mea- 
sured pressure distribution for a high Reynolds number with the 
inviscid pressure distribution on the blade. In addition to the up- 
stream bleed flaps, a tailboard on the pressure side was used to set 
the pressure gradient. The compiled data include pressure surveys, 
mean and fluctuating velocity profiles, intermittency profiles, and 
turbulent shear stress profiles. It was observed that the effect of 
high Reynolds number or high freestream turbulence level was to 
move transition upstream. Transition started in the shear layer 
over the separation bubble and led to rapid boundary layer reat- 
tachment. At the lowest Re case, transition did not take place 
before the trailing edge and the boundary layer did not reattach. 
The beginning of transition corresponded to the beginning of a 
significant rise in the turbulent shear stress. These experimental 
results provide detailed documentation of the boundary layer and 
extend the existing database to lower Reynolds numbers. The 
cases used for comparisons with computations in this paper are 
listed in Table 1 along with the type of data used for comparisons. 


3 Intermittency Transport Model 

In this section, the transport model for intermittency is pre- 
sented. The model combines the transport equation models of 
Steelant and Dick [9] and Cho and Chung [10]. Details of the 
development and implementation of the transport model are given 
in Suzen and Huang [11,16,17], and in Suzen et al. [18]. 

The model equation is given by 
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The distributed breakdown function, /(s) has the form 
as' A + bs' 3 + cs' 2 + ds' +e 


f(s) 


gs' 3 + h 


( 1 ) 


( 2 ) 


where s' =s-s„ and s is the distance along the streamline coordi- 
nate, and s, is the transition location. The coefficients are 


a = 50\ — b = - 0.4906 
U 


c 


0.2041 



d = 0.0 


e = 0.044441 


ncr 

U 


- 1.5 


h= We 

The shear stresses are defined as 


T a = H 


du t du^ 
dXj dx t 


J 


g = 50 


(3) 

2 du k 

2 

(4) 

3 dej . 

- -pkSij 


The blending function F is constructed using a nondimensional 
parameter k/Wv, where k is the turbulent kinetic energy and W is 
the magnitude of the vorticity. The blending function has the form 


F = tanh 4 


kJWv 


L2oo(i - y )i )°- 3 j 

The model constants used in Eq. (1) are 

Cyi = (Ty, = 1 .0 C 0 =1.0 C) = 1.6 


(5) 


C 2 = 0.16 C 3 = 0.15 

The intermittency is incorporated into the computations simply by 
multiplying the eddy viscosity obtained from a turbulence model, 
m, by the intermittency factor, y. Simon and Stephens [4] showed 
that, by combining the two sets of conditioned Navier-Stokes 
equations and making the assumption that the Reynolds stresses in 
the nonturbulent part are negligible, the intermittency can be in- 
corporated into the computations by using the eddy viscosity, fi* 
which is obtained by multiplying the eddy viscosity from a turbu- 
lence model, p, t , with the intermittency factor, y. That is 

Ft = yFt (6) 

is used in the mean flow equations. It must be noted that y does 
not appear in the generation term of the turbulent kinetic energy 
equations. 

Computations of the experiments are performed using a re- 
cently developed multiblock Navier-Stokes solver, called ghost. 
The code was developed at the University of Kentucky, by Huang, 
and is a pressure-based code based on the SIMPLE algorithm with 
second-order accuracy in both time and space. Advection terms 
are approximated by a QUICK scheme and central differencing is 
used for the viscous terms. The “Rhie and Chow” momentum 
interpolation method [25] is employed to avoid checkerboard os- 
cillations usually associated with the nonstaggered grid arrange- 
ment. This code is capable of handling complex geometries, mov- 
ing, and overset grids and includes multiprocessor computation 
capability using message passing interface (MPI). Since multiple 


processors are used during the computations, it is more efficient to 
divide the computational domain into several smaller pieces with 
very fine grids and distribute the zones to processors with the 
consideration of load balancing. This code has been used exten- 
sively in a recent turbulence model validation effort (Hsu et al. 
[26]) and computations of unsteady wake/blade interaction (Suzen 
and Huang [27]) conducted at the University of Kentucky. 

The multiblock grid systems used in the computations are ob- 
tained by conducting a series of grid refinement studies in order to 
ensure that the details of the flow field are captured accurately and 
the results are grid independent. All grid systems have first y + less 
than 0.5 near solid walls. 

In using this intermittency approach, the turbulence model se- 
lected to obtain p. t must produce fully turbulent features before 
transition location in order to allow the intermittency to have full 
control of the transitional behavior. Menter’s [28] SST model sat- 
isfies this requirement. It produces almost fully turbulent flow in 
the leading edge of the boundary layer and therefore is used as a 
baseline model to compute p. t and other turbulent quantities in the 
computations [18], 

The value of ncr used in evaluating the constants given by Eq. 
(3) is provided by the following correlation for zero-pressure gra- 
dient flows [18] 

n<r = (nv 2 /U 3 )cr= 1.8 X 10- n 7n 7/4 (7) 

When flows are subject to pressure gradients, the following cor- 
relation is used 

her _ f M^ 1 - exp (°- 75 * io 6 x,7'ir 0 - 7 )] , *,<() 

(ncrW" [i 0 -322< 5985 , K,>0 (8) 

with M defined as 

M = (8507m- 3 - 1007V 0 - 5 + 120) 

where (ner) ZPG is the value for flow at zero pressure gradient and 
can be obtained from Eq. (7), and K,=(v/U 2 )(dU/dx) t is the flow 
acceleration parameter. The favorable pressure gradient part of the 
above correlation (for K,> 0) is from Steelant and Dick [9], The 
portion of the correlation for adverse pressure gradient flows for 
K,< 0 is formulated using the transition data of Gostelow et al. [6] 
and Simon et al. [14] (Suzen et al. [18]). 


1 

0.5 

U 0 

>. 

- 0.5 
-1 

Fig. 6 Multiblock grid used for computations of experiments 
of Lake et ai. [3,22] and FSTI=0.08% experiments of Huang et 
al. [23] 



i . i i i 

- 0.5 0 0.5 1 1.5 

*/C x 


Journal of Turbomachinery 


JULY 2007, Vol. 129 / 531 


Downloaded 10 Dec 2008 to 128.156.10.80. Redistribution subject to ASME license or copyright; see http://www.asme.org/terms/Terms_Use.cfm 







Fig. 7 Comparison of computed pressure coefficient with experiments of Lake et al. [3,22] 


The current approach uses the intermittency transport model to 
obtain the intermittency distribution for the transitional flows, 
while the onset of transition is defined by correlations. 

The onset of attached flow transition is determined by the fol- 
lowing correlation in terms of turbulence intensity, Tu, and the 
acceleration parameter, K„ 

Re ft = (120 + 1507’«- 2/3 )coth[4(0.3 -K,X 10 5 )] (9) 

where K t was chosen as the maximum absolute value of that pa- 
rameter in the downstream deceleration region [18]. This correla- 
tion maintains the good features of Abu-Ghannam and Shaw [29] 
correlation in the adverse pressure gradient region, and in addition 
reflects the fact that the flow becomes less likely to have transition 
when subject to favorable pressure gradients by rapidly rising as 
K, becomes positive. 

In order to determine the onset of separated flow transition Re s( 
is expressed in terms of the turbulence intensity (Tu) and the 
momentum thickness Reynolds number at the point of separation 
(Re ft ) in the form [19] 

Re s( = 874Re^ 71 exp[- 0.47 m] (10) 

This correlation provides a better representation of the experimen- 
tal data than Davis et al. [30] correlation and is used to predict 
onset of separated flow transition in the present computations. 
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4 Results and Discussion 

4.1 Simulations of Experiments of Lake et al. [3,22]. The 

intermittency model is applied to predict the Pack B blade experi- 
ments of Lake et al. [3,22], In the computations, flows at Rey- 
nolds numbers of 86,000 and 172,000 based on inlet velocity and 
axial chord with ffeestream intensities of 1% and 4% were 
investigated. 

The computations were performed using the grid system shown 
in Fig. 6 consisting of five zones obtained as a result of a grid 
refinement study. In the grid refinement study computations were 
performed on a series of successively finer grids and the variations 
in the results were observed. The grid shown in Fig. 6 was chosen 
to be adequate for obtaining grid-independent solutions for all 
cases. The four zones on which the blade grid is superposed each 
have 125 X 225 grid points and the O-type grid around the blade 
has 401 X 101 points with first y + less than 0.5. 

The comparisons of computed and experimental pressure coef- 
ficient distributions are shown in Figs. l(a)-l(d). In these figures, 
the experimental distributions correspond to the measurements 
made on test blades 4 and 6. 

The computed results compare well with the experiments for 
high turbulence intensity, FSTI=4%, cases shown in Figs. 7(a) 
and 7(c). However, for FSTI=1% cases shown in Figs. 1(b) and 
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50000 100000 150000 200000 

Re (U ln C x /v ) 


Fig. 8 Comparison of computed total pressure loss coeffi- 
cients with experiments of Lake et al. [3,22] and Huang et al. 
[23] 


1(d), the extent of the separation bubbles is underpredicted in the 
computations. For example, for Re =86,000, FSTI=1%, shown in 
Fig. 1(d), the flow reattaches earlier in computations than it does 
in the experiment, as can be observed from the difference in the 
pressure coefficient distributions between x/C x =&0 to 85%. 

The comparison of computed total pressure loss coefficients 
with experiments is shown in Fig. 8. For the Re= 86,000 case, the 
computed loss coefficient is in good agreement with the experi- 
ments for both FST1 levels. However, for the Re= 172,000 case 
the computations underpredicted the loss coefficient compared to 
experiments for both FSTI=1% and FSTI=4%. From Fig. 8 it is 
evident that the cascade losses decrease as the Reynolds number 
increases. This reduction in cascade losses with increasing Rey- 
nolds number is due to the decrease in size of the separated flow 
region on the suction side of the blades. 

The onset of separation locations, reattachment locations, and 
onset of transition locations on the suction surface are summa- 
rized in Table 2 for these cases, along with the corresponding 
values from experiments. In the experiments, the onset of transi- 
tion locations and the reattachment locations are not reported. The 
experimental onset of separation and reattachment points are ex- 
tracted from the experimental pressure coefficient data. The onset 
of separation is taken to be the axial location where the plateau in 
the pressure coefficient distribution of the suction side begins, and 
the reattachment point is taken to be the axial location after the 
sharp change in C p following the plateau. This procedure may 
lead to an error of approximately ±1.5% of axial chord in the 
estimated onset locations. 

The onset of separation, reattachment, and onset of transition 
locations are plotted against Reynolds number in Figs. 9(a) and 
9(b) for FSTI=4% and 1%, respectively. The uncertainty in the 
estimated experimental values is indicated by error bars in the 
figures. For the high turbulence intensity case, computation pre- 
dicts onset of separation and reattachment slightly upstream of the 
experiment. For the low FSTI case shown in Fig. 9(b), the sepa- 
ration zone is predicted smaller than the experiments. The onset of 
transition is predicted over the separated flow region in the shear 



50000 100000 150000 200000 

Re(U, n C„/v) 

Fig. 9 Comparison of separation, reattachment, and transition 
locations for experiments of Lake et al. [3,22] 


layer. From comparison of these figures it is evident that, with 
decreasing freestream turbulence intensity, the separation zone be- 
comes larger, and for a given FSTI condition, the separated flow 
region gets smaller with increasing Reynolds number. 

4.2 Simulations of Experiments of Huang et al. [23]. In this 
set of experiments, first the cases with no grid in tunnel corre- 
sponding to FSTI =0.08% are computed. In these computations, 
the same grid system used for the computations of experiments of 
Lake et al. [3,22] shown in Fig. 6 is used. 

The comparisons of the computed and the experimental pres- 
sure coefficients are shown in Figs. 10(a)-10(e) for Re= 100,000, 
75,000, 50,000, 25,000, and 10,000 based on inlet velocity and 
axial chord. The agreement between the experiments and compu- 
tations is very good for all cases. 

The computed total pressure loss coefficients are compared to 
the available data for Re=25,000 and 50,000 in Fig. 8. The loss 
coefficients predicted in the computations are 2% to 3% higher 
compared to the experiments for both Reynolds numbers. 

The onset of separation, transition, and reattachment locations 
are tabulated in Table 3 for all cases and plotted against Reynolds 
number in Figs. ll(a)-ll(c) for FSTI=0.08%, 1.6%, and 2.85%, 


Table 2 Separation, reattachment, and transition locations for cases of Lake et al. [3,22] 
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Fig. 10 Comparison of computed pressure coefficients with experiments of Huang et al. [23] for FSTI=0.08% cases 


respectively. 

Computed velocity profiles at seven axial stations along the 
suction surface of the blade are compared to the experiments for 
Re= 100,000, 75,000, and 50,000 in Figs. 12-14, respectively. 

For the Re= 100,000 case, the computed velocity profiles com- 
pare very well with the experiment as shown in Figs. 12(a)-12(g). 
At the first three measurement stations, flow is laminar and at- 
tached as shown in Figs. 12(a)-12(c). Flow separation takes place 
at x/C x = 0.725 and the separated flow region is visible in Figs. 
12(d) and 12(e), corresponding to axial locations of x/C x = 0.75 
and 0.80. The flow transition and reattachment takes place around 
x/C x =0.84 in the computation. Reattachment location is earlier 
than the experiment which takes place at x/C x =0.875. In Fig. 
12 (J) corresponding to axial station of x/C x =0.85 the computed 
flow field has already attached, although the experimental profile 
indicates a very small separation zone close to wall. At xlC x 
=0.9 the flow is completely attached as shown in Fig. 12(g). 

When the Reynolds number is reduced to 75,000, the size of the 
separation bubble increases as can be observed from the compari- 
son of the velocity profiles shown in Figs. 13(a)-13(g). At this 
Reynolds number the flow separates around x/C x ~ 0.72 and reat- 
taches around x/C x ~0.87. The transition onset location is pre- 


dicted at x/C x = 0.854. The size of the separation bubble is larger 
than the Re= 100,000 case from comparison of Figs. 13(d)-13(/) 
and 12(t/)-12(/). 

Next, the Reynolds number is reduced to 50,000 and the com- 
parison of computed and experimental velocity profiles is shown 
in Figs. 14(a)-14(g). For this case the separation bubble is much 
larger from the previous cases and extends until x/C x ~ 0.975 in 
the experiment and x/C x ~0.93 in the computations, as can be 
seen in Figs. 14(d)-14(g). Computations predicted the transition 
onset location at x/C x = 0.89. In the computations, the onset of 
separation is predicted well in agreement with experiment; how- 
ever, the reattachment point is earlier, making the size of the sepa- 
ration bubble smaller when compared to experiment. This is evi- 
dent from the comparison of velocity profiles at the last two 
stations shown in Figs. 14 (/) and 14(g). 

The onset of separation and reattachment points for FSTI 
=0.08% cases is predicted upstream of the experiments as shown 
in Fig. 1 1 (a). 

Next, the high FSTI cases are computed using the six zone 
multiblock grid system shown in Fig. 15. The computational do- 
main is extended upstream of the blade in order to specify the 
correct turbulence intensity at the inlet and to match the decay of 


Table 3 Separation, reattachment, and transition locations for cases of Huang et al. [23] 


Re 

GW*) 

FSTI 

(%) 

x,IC x 

(Computation) 

x,IC x 

(Experiment) 

x r /C x 

(Computation) 

x r IC x 

(Experiment) 

x trlCx 

(Computation) 

10,000 

0.08 

0.661 

0.725 




25,000 

0.08 

0.656 

0.725 

0.980 


0.936 

50,000 

0.08 

0.714 

0.725 

0.925 

0.975 

0.890 

75,000 

0.08 

0.718 

0.725 

0.860 

0.870 

0.854 

100,000 

0.08 

0.725 

0.725 

0.840 

0.875 

0.840 

50,000 

1.6 

0.722 

0.728 

0.900 

0.900 

0.854 

75,000 

1.6 

0.728 


0.867 

0.875 

0.834 

100,000 

1.6 

0.732 


0.860 

0.877 

0.821 

50,000 

2.85 

0.728 

0.722 

0.887 

0.900 

0.837 

75,000 

2.85 

0.732 

0.729 

0.840 

0.870 

0.816 

100,000 

2.85 

0.735 

0.734 

0.842 

0.850 

0.806 
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Fig. 11 Comparison of separation, reattachment, and transi- 
tion locations for experiments of Huang et al. [23] 


turbulence that reaches the blade. The matched computed and ex- 
perimental turbulence decays are shown in Figs. 3 and 4 for grid 
0 and grid 3, respectively. The cases considered have the grids 
placed 0.762 m (30 in) upstream of the blade, corresponding to 
turbulence intensities of 2.85% and 1.6% for grid 0 and grid 3, 
respectively. 

The comparison of the computed and the experimental pressure 
coefficient distributions for Re=50,000, 75,000, and 100,000 for 
FSTI=2.85% cases is shown in Fig. 16. The agreement is very 
good between computations and experiments. 

Comparisons of computed velocity profiles with the experi- 
ments for Re= 100,000 are given in Figs. 17(a)-17(g). In this 
case, the flow separates around x/C x = 0.74 and reattaches at 
x/C x =0.85. The onset of transition is predicted at x/C x =0.806. 

Journal of Turbomachinery 


The computed size and extent of the separation bubble is in good 
agreement with the experiment as tabulated in Table 3 and as can 
be seen in Figs. 11(c) and 16(d)-16(j). 

For the lower Reynolds number of 75,000, computed velocity 
profiles are compared with the experiments in Figs. 18(a)-18(g). 
The agreement between experiment and computation is good prior 
to the reattachment as shown in Figs. 18(a)-18(e). There is a 
discrepancy in the reattachment region. The flow separation takes 
place around x/C x ~0.73 and reattaches at x/C x ~0.87 according 
to the experiment, whereas computation predicts reattachment ear- 
lier at around x/C x ~0.84 with the onset of transition predicted at 
x/C x =0.816. The difference in reattachment points is evident in 
the comparison of the computed and experimental velocity pro- 
files shown in Fig. 18(/). At this station the experimental profile 
indicates separated flow and the computed profile shows an al- 
ready attached flow. 

The next case considered has the same FSTI=2.85% but with 
Reynolds number being reduced to 50,000. The comparison of 
velocity profiles is shown in Figs. 19(a)-! 9(g). The computations 
agree well with the experiment, and the size and extent of the 
separation bubble are well predicted as can be seen from Fig. 
11(c). The onset of separation is around x/C x ~0.72 and the flow 
reattaches around x/C x ~0.9, with transition onset at x/C x 
=0.837. 

In Fig. 20, computed and experimental pressure coefficient dis- 
tributions for grid 3 case which correspond to FSTI=1.6% are 
compared for Re=50,000, 75,000, and 100,000. Again, very good 
agreement between computations and experiments is obtained. 
The onset of separation and reattachment locations shown in Fig. 
11(h) compares well with the experiments. 

Overall, Figs. ll(a)-ll(c) indicate that, as FSTl increases, the 
separated flow region decreases, and at a given FSTI, increasing 
Reynolds number has the same effect on the separated flow re- 
gion. 

4.3 Simulations of Pack B Experiments of Volino [24], In 

computation of experiments of Volino [24] the flow field is mod- 
eled with the 31 -zone multiblock grid shown in Fig. 21 obtained 
as a result of a series of grid refinement studies. The bleed flaps 
below the lower blade and above the upper blade are defined by 
fitting third-order polynomials through the available points ob- 
tained from experimental setup; these curves are used as the flap 
shapes in generating the computational grid. Initial computations 
indicated that the shape of the bleed flaps and the orientation of 
the tailboard behind the upper blade greatly affect the computed 
results, especially the onset of separation and reattachment points 
on the lower blade’s suction surface. In order to select the most 
accurate orientation for the tailboard and the shape of the bleed 
flaps, several test computations were performed for the case with 
Re=41,162 and FSTI=0.5% using different tailboard orientations 
and bleed flap shapes. In these computations the main goal was to 
match the experimental velocity profiles in the laminar flow part 
and to capture the correct onset point of separation. Once an ac- 
ceptable geometry is obtained, the final bleed flap shapes and 
tailboard orientation are used for computation of all other Rey- 
nolds number cases. 

Computed pressure coefficient distributions are compared to 
experiments in Figs. 22(a)-22(d) for Re=82,324, 41,162, 20,581, 
and 10,291, and the separation onset, reattachment, and transition 
onset information is summarized in Table 4. The C p comparison 
for Re = 82 , 324 shown in Fig. 22(a) indicates that the computation 
predicts early reattachment of the flow; in the recovery region 
following reattachment the pressure coefficient distribution is 
overpredicted. 

The computed pressure coefficient distributions for the lower 
Reynolds number cases shown in Figs. 22(c) and 22(d) compare 
well with experiments. For the Re=41,162 case shown in Fig. 
22(h), the onset of separation and reattachment locations matches 
the experiment as given in Table 4; however, in the recovery re- 
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Fig. 12 Comparison of computed velocity profiles with experiments of Huang et al. [23], Re=100,000, FSTI=0.08% case 



Fig. 13 Comparison of computed velocity profiles with experiments of Huang et al. [23], Re=75,000, FSTI=0.08% case 



Fig. 14 Comparison of computed velocity profiles with experiments of Huang et al. [23], Re=50,000, FSTI=0.08% case 


Ih 



- 2-101 


x/C x 

Fig. 15 Grid used for computation of experiments of Huang et 
al. [23] with FSTI=1.6% and 2.85% 


gion the pressure coefficient distribution is overpredicted. 

Computed velocity profiles are compared to experiment at 1 1 
stations along the suction surface of the blade in Figs. 23(a)-23(k) 
for Re=82,324 and FSTI=0.5%. The results compare well with 
the experiment up to x/C x = 0.732 shown in Figs. 23(a)-23(g). 
After this station flow separation takes place. Separation onset and 
reattachment are slightly earlier in the computations compared to 
experiment as given in Table 4. This also can be observed from 
the velocity profiles at stations xl C x =0.798 to 0.912 shown in 
Figs. 23(h)-23(j). Overall computations compare well with the 
experimental measurements. 

Next the Reynolds number is reduced to 41,162 and the com- 
puted and experimental velocity profiles are compared in Figs. 
24(a)-24(k). The computed profiles agree well with experiments 
except at jc/C x = 0.912 shown in Fig. 2A(j). At this station the 
computation indicates a smaller separated flow region close to 
reattachment in contrast to the experiment. However, the flow 
reattaches around x/C x = 0.95 both in computation and experi- 
ment, and in the next measurement station the agreement is well. 

The next case considered has a Reynolds number of 20,581. 
Computed velocity profiles are shown along with the experimental 
data at 11 axial stations in Figs. 25(a)-25(k). In this case flow 
separates around xl C x ~ 0.76 and does not reattach in experiment; 
however, computations indicated reattachment at xl C X ~Q.9%. 
This discrepancy is evident from the comparison of velocity pro- 
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Fig. 16 Comparison of computed pressure coefficients with experiments of Huang et al. [23] for FSTI=2.85% cases 
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Fig. 17 Comparison of computed velocity profiles with experiments of Huang et al. [23], Re=100,000, FSTI=2.85% case 
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Fig. 18 Comparison of computed velocity profiles with experiments of Huang et al. [23], Re=75,000, FSTI=2.85% case 
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Fig. 19 Comparison of computed velocity profiles with experiments of Huang et al. [23], Re=50,000, FSTI=2.85% case 
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Fig. 20 Comparison of computed pressure coefficients with experiments of Huang et al. [23] for FSTI=1.6% cases 
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Fig. 21 Thirty-one zone multiblock grid used for computation 
of experiments of Volino [24] 




files at the last two measurement stations shown in Figs. 25 (j) and 
25 (k). The computation indicates a smaller separated region in 
these stations and finally reattaches very close to the trailing edge. 
Onset of transition was predicted at xlC x = 0.978. 

The final case in this set of experiments is the one with Re 
= 10,291. The computed velocity profiles compare very well with 
the experimental data as shown in Figs. 26(a)-26(k). In this case 
the flow separates around x/C x ~0J6 and does not reattach. The 
flow is completely laminar; transition was not predicted on the 
blade. 


5 Concluding Remarks 

A transport equation for the intermittency factor is employed to 
predict three sets of recent low-pressure turbine experiments on 
the Pack B blade. The intermittent behavior of the transitional 
flows is taken into account by modifying the eddy viscosity with 
the intermittency factor. Comparisons of the computed and experi- 
mental data are made and overall good agreement with the experi- 
mental data is obtained. The predicting capabilities of the current 
intermittency approach and the intermittency transport model in 
prediction of transitional flows under a wide range of low- 
pressure turbine conditions is demonstrated. 




Fig. 22 Comparison of computed pressure coefficient distributions with experiments of Volino [24], FSTI=0.5% 
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Table 4 Separation, reattachment, and transition locations for cases of Voiino [24] 


Re 

(U in C x /v) 

FSTI 

(%) 

xj C x 

(Computation) 

x s IC x 

(Experiment) 

x,/C x 

(Computation) 

x r /C x 

(Experiment) 

x tr! Cx 

(Computation) 

10,291 

0.5 

0.760 

0.750 




20,581 

0.5 

0.765 

0.760 

0.980 


0.978 

41,162 

0.5 

0.760 

0.770 

0.950 

0.950 

0.840 

82,324 

0.5 

0.757 

0.767 

0.890 

0.900 

0.857 


x / C x = 0 093 |(x/C x = 0.J97 | [_ x/ C x = 0.3j>14 | [ x / C x = 0.4(35 ||_x/C x = 0.5E 

Computation 

o Experiment 


(f) 

x/C =0. 


: 

(a) 

x / C K = 0 

093 

Compulation 

o Experiment 




(b) 

x/C =0.197 


■(c) 

x/C =0.3-14 


x/C =0.73 


(h) 

x/C =0.79 


Fig. 24 Comparison of computed velocity profiles with experiments of Voiino [24], Re=41,162, FSTI=0.5% 
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Nomenclature 

C p = pressure coefficient, 2(P-P 00 )/(p 00 U? n ) 

C x = axial chord 

FSTI = freestream turbulence intensity (%) 

K, = flow acceleration parameter (v/U*)(dU e /ds) 
k = turbulent kinetic energy 
L x = axial chord 

N = nondimensional spot breakdown rate param- 
eter, n<T(? t / v 
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n = spot generation rate 
P = static pressure 
P tot ai = total pressure 
Re = Reynolds number 
R e J( = ( s t -s s )U e /v 
Reg, = 8,U e /v 

s = streamwise distance along suction surface 
Tu = turbulence intensity (%), u' I U 
U = boundary layer streamwise velocity 
U e = local freestream velocity 
U m = inlet freestream velocity 
u T = friction velocity 
W = magnitude of vorticity 
y n = distance normal to the wall 

y + = y„ujv 

y = intermittency factor 
8 = momentum thickness 
\g = pressure gradient parameter (8 1 / v)(dU/ds) 

/ x i = molecular viscosity 
fju, = eddy viscosity 

v = flip 

v, = fi t /p 

o) = total pressure loss coefficient, 

2 ( P total inlet -^totalJ I 

p = density 

a = spot propagation parameter 

Subscripts 

e = freestream 
s = onset of separation 
t = onset of transition 
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